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The consequences of parity symmetry for higher-order statistics 
of cosmic shear and other polar fields 
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Abstract. We consider the parity transformation and the consequences of parity invariance on the n-point cor- 
relation function of the cosmic shear caused by gravitational lensing by the large-scale structure, or any other 
polar (or 'spin-2') field. The decomposition of the shear field into E- and B-modes then yields the result that 
any correlation function which contains an odd number of B-mode shear components vanishes for parity-invariant 
random shear fields. In particular, this result implies that the expectation value of the third-order cross aperture 
statistics, (M*), vanishes for parity-invariant shear fields. Therefore, a significant detection of a non-zero value 

of (M x ) in a cosmic shear survey does not indicate the presence of B-modes, but of an underestimate of the 
statistical uncertainty or cosmic variance, or a remaining systematic effect in the shear measurement. We argue 
that the parity invariance provides a very specific diagnostic of systematic effects in shear data. Our results apply 
as well to the linear polarization of the cosmic microwave background. 
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1. Introduction 

Cosmic shear, resulting from the gravitational deflection 
of light bundles from distant sources by the tidal grav- 
itational field of the large-scale structure, has been rec- 
ognized as an important tool for observational cosmology 
(e.g., Blandford et al. 1991; Miralda-Escude 1991; Kaiser 
1992; see Mellier 1999 and Bartelmann & Schneider 2001 
for recent reviews). It allows the direct investigation of 
the (dark) matter distribution up to redshifts ~ 1 , without 
making assumptions about the relation between dark mat- 
ter and luminous tracers, like galaxies. Therefore, since the 
first detections of cosmic shear (Bacon et al. 2000; Kaiser 
et al. 2000; van Waerbeke et al. 2000; Wittman et al. 2000) 
intensive observational and theoretical research has been 
performed. 

Most of the work up to now has concentrated 
on second-order statistics of the shear field. However, 
Bernardeau et al. (2002) have detected a significant third- 
order shear signal in the Virmos-Descart cosmic shear sur- 
vey. It was realized before that the third-order cosmic 
shear statistics would provide an invaluable tool for cos- 
mological parameter estimates; in particular, a third-order 
shear measurement breaks the (near) degeneracy between 
the density parameter Q m and the power spectrum nor- 
malization erg (Bernardeau et al. 1997; Jain & Seljak 1997; 
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van Waerbeke et al. 1999). These theoretical studies con- 
centrated on the third-order statistics of the convergence 
field, from which the shear field is derived. However, as 
the shear is the observable (or more precisely, the observ- 
able galaxy ellipticities provide an unbiased estimate of 
the shear), not the convergence, recent interest has fo- 
cused on determining the third-order statistical properties 
of the shear itself. 

Since the shear is a two-component quantity, the three- 
point correlation function (3PCF) has eight independent 
components and depends on three arguments. Bernardeau 
et al. (2003) considered a specific combination of com- 
ponents of the 3PCF, calculated its properties, both an- 
alytically and numerically, and applied this to measure 
the 3PCF on the Virmos-Descart data (Bernardeau et al. 
2002). With the large number of components of the 3PCF, 
it is not clear a priori which of them carry most of the in- 
formation about the underlying matter field. In contrast 
to the two-point correlation function (2PCF), which has 
only two non- vanishing components (the other two vanish 
by parity invariance) , none of the eight components of the 
3PCF vanishes for a general set of three points. Only for 
specific triangles (where two or three sides are equal) does 
parity imply that some components vanish (Schneider & 
Lombardi 2003; hereafter SL03). This has later been veri- 
fied using ray-tracing simulations through a cosmological 
matter distribution, as well as with analytic estimates of 
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the 3PCF from the halo model (Takada & Jain 2003a; 
Zaldarriaga & Scoccimarro 2003). As it turned out, inves- 
tigating the statistical properties of the 3PCF of a polar 
field is far from trivial, and some confusion about symme- 
try properties has arisen. 

The purpose of this paper is to clarify some of the is- 
sues related to parity transformations of the shear field. In 
Sect. 2 we recall the definition of the shear 3PCF and its 
behavior under parity transformations. In Sect. 3 the de- 
composition of a general shear field into E- and B-modes 
(Crittenden et al. 2002; Schneider et al. 2002) is intro- 
duced. The central result of this paper, derived in Sect. 4, 
is as follows: For a pure B-mode shear field whose statis- 
tical properties are invariant under parity transformation 
(hereafter: parity symmetric), the 3PCF vanishes identi- 
cally. More generally, the n-point shear correlation func- 
tion containing m B-mode shear components and (n — m) 
E-mode shear components vanishes identically for m odd. 
This result will be shown both by an explicit calculation, 
as well as by a simple heuristic, though accurate, consid- 
eration (at the beginning of Sect. 4). 

2. Parity transformation of the 3PCF 

2.1. Definition of the 3PCF 

We consider a shear field (or more generally, any polar 
field) that is statistically homogeneous and isotropic. This 
implies that the 3PCF depends on three arguments only, 
which define the shape and size of a triangle. Let Xj, i = 
{1, 2, 3}, be three points defining a triangle, and let xi = 
X3 — X2, X2 = Xi — X3, and X3 = X2 — Xi be the 
side vectors. Furthermore, let 7 M (Xj), fi = {1,2}, be the 
Cartesian components of the shear at the point Xj. We 
note that the Cartesian components of the shear can also 
be considered as components of a symmetric trace-free 
matrix, 



r = 



71 72 

72 -71 



(i) 



and its behaviour under linear coordinate transformations 
x — > x' = Tx is 

T' = T T TT. (2) 
Then, the Cartesian components of the 3PCF are defined 



7 aii ,a(xi,x 2 ,x 3 ) = (7^X1) 7^X2) 7a(X 3 )} 



(3) 



where it was used that the shear field is a homogeneous 
random field, so that the 3PCF is invariant under transla- 
tions; therefore, 7 jUI /a does depend only on the separation 
vectors x;. In order to employ the assumed rotational in- 
variance of the shear field, one needs to measure the shear 
components relative to a direction that rotates with the 
triangle. For example, one can define the center of mass 
o = (Xi + X2 + Xs)/3 of the triangle and define the shear 
components at Xj relative to the direction vector Xj — o 
between the point Xj and the center o. If Q denotes the 



polar angle of the direction vector Xj — o, we define the 
tangential and cross component of the shear at Xj relative 
to the center o by 

7(Xj,Ci) = 7 t(X l ,0) + i7x(X l ,C i ) := - 7 (X0e- 2i<: « 

= -[ 7 i(X l )+i72(X J )]e- 2i <* , (4) 

as can be seen from Q when setting T as a rotation ma- 
trix. We then define the components of the 3PCF relative 
to the center of mass as 

jj^X ( x U^2,x 3 ) = (7m( x i)7^(X 2 )7a(X 3 )) , (5) 

where now /j,,v, A = {t, x}. It should be noted that the 
3PCF defined in this way is invariant under rotations of 
the triangle, and thus depends only on the lengths Xj = 
|xj| of the side vectors, provided that we fix the orientation 
such that Xi x X2 = X2 x X3 = X3 x xi > 0. 

There is nothing special about the center of mass; in- 
stead, different centers of the triangle (or other points at- 
tached to a triangle) can be chosen in the definition of the 
tangential and cross components of the shear, and thus 
in defining the 3PCF; SL03 considered some of the most 
obvious choices (orthocenter, center of incircle, center of 
circumcircle) , and derived the relations between the 3PCF 
defined with respect to different centers. 

2.2. Parity transformation 

SL03 have derived the behavior of the 3PCF under par- 
ity transformations. Consider a fixed triangle, and its 
mirror symmetric sister obtained by flipping the origi- 
nal one along the line perpendicular to the side vector X3 
going through X3. This transformation has two effects: 
first, the points Xi and X2 are interchanged, so that the 
orientation of the triangle is now negative, in the sense 
that xi x X2 < 0. Second, the tangential components of 
the shears are unaffected, whereas the cross components 
change sign. Together, this implies for the behavior under 
parity transformations that 

nT;™'^,^,^), (6) 



(ccn) t \ 



where n is the parity; it is +1 if all of the indices of 7 are 
t's, or two x's occur, otherwise it is —1. Those compoents 
of the shear 3PCF for which n = +1 are called even, 
the others are the odd components. Note the interchange 
of the arguments in ©. Furthermore, 10 applies equally 
well to the 3PCF of the shear if measured with respect to 
a different reference point attached to the triangle. 

We can restate this result in the notation of Takada 
& Jain (2003a, b). Instead of specifying a triangle by its 
side lengths, one can also characterize it by two sides 
and one angle. Let ip be the angle at the vertex X3 be- 
tween the directions to Xi and X2, where ip £ [0, it] 
if Xi x X2 > 0, and ip G [tt, 2ir] otherwise. Then, 
1 ft vx{xi,x 2 ,ip) = Jpv\(xi,X2,x 3 ) for ip G [0,7r], and 
7nv\(xi, x 2 ,ip) = Ju li x(x 2 ,Xi,X3) for ip e [tt, 27r] in the 
new notation, and the parity transformation J(|J) becomes 

P hnux{xi,x 2 ,ip)\ = Ujfj,„\(xi,x 2 ,2iT - ip) . (7) 
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3. E- and B-modes of the shear 

Provided the shear field is due to gravitational lensing by a 
(geometrically-thin) matter distribution, it can be derived 
from a deflection potential ip, according to 



7 = 7i + 172 



2 (^,11 _ ^,22) + 1^,12 



(8) 



where indices separated by a comma denote partial deriva- 
tives with respect to the Cartesian coordinates. A shear 
field with the property JHJ is called an E-mode field. A 
general shear field will not satisfy JSJl, i.e., will not be 
derivable from a single scalar potential. However, a gen- 
eral shear field can be written as (see Crittenden et al. 
2002; Schneider et al. 2002; see also Bunn et al. 2003 for 
an E/B-mode decomposition of the CMB polarization) 
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f 



f 



V^12 + o (V£l " ^22) 



2 tyjn - ^22) - V> B i2 

= {if + 7?) + i (7 2 E + 7 2 B ) ■ 0) 

with the two potentials ip E - B , which are conveniently com- 
bined into a complex field 



One also defines the Laplacians of the potentials, 



2„/,E,B 



2k 



E,B 



(10) 



(11) 



in lensing, k e is the convergence (or dimensionless sur- 
face mass density) of the deflector. In the case of a matter 
distribution extending between sources at high redshift 
and us, as for lensing by the large-scale mass distribu- 
tion in the Universe, the geometrically-thin lens theory 
no longer applies; however, both analytical considerations 
(e.g., Bernardeau et al. 1997; Schneider et al. 1998; van 
Waerbeke et al. 1999), as well as numerical investigations 
(e.g., Jain et al. 2000) have shown that the cosmic shear 
field is very close to an E-mode field; the B-mode contribu- 
tions coming from multiple deflections are suppressed by a 
large factor compared to the E-mode shear. Furthermore, 
a small B-mode contribution derives from the spatial clus- 
tering of source galaxies (Schneider et al. 2002), but again 
is very much smaller than the E-mode, except at very 
small angular scales. However, if the shear is due not only 
to lensing, but also to, e.g., intrinsic alignments of galax- 
ies (e.g., Heavens et al. 2000, Crittenden et al. 2001; Croft 
& Metzler 2000; Catelan et al. 2000), a non- vanishing B- 
mode can be expected. Significant B-mode contributions 
to the cosmic shear signal have been measured (e.g., van 
Waerbeke et al. 2001a, 2002; Hoekstra et al. 2002; Jarvis 
et al. 2003). 

It should be noted that the parity transformation be- 
havior © or (JZJ) are valid for a general shear field; in 
the derivation, no assumption on the E/B-mode charac- 
ter of the shear field has been made. Also, we have not 
yet used the fact that the statistical properties of the 
shear field are invariant under parity transformation (i.e., 
parity-symmetric). Takada & Jain (2003a, b) have argued 



that the transformation laws ©, 10 only apply for E- 
modes, and that the signs are changed for a pure B-mode 
field. Their argument shall be reproduced here: 

Consider a pure E-mode field, for which applies. 
This can be turned into a B-mode field by rotating all 
shears by 45 degrees; in this operation, 7t — > 7 X , and 
7x — > — 7t- Hence, what was 7ttt for the E- field becomes 
7x x x for the B-mode field, implying that 7x x x does not 
change sign under a parity transformation. This argument, 
however, is incomplete: the parity-reversed B-mode shear 
field is obtained from the parity-reversed E-mode shear 
field by rotating the shear by —45 degrees, and this guar- 
antees that remains valid even for a pure B-mode field. 
As mentioned above, the derivation of JH| and makes 
no assumption about the character of the shear field, but 
is derived purely from geometrical considerations (liter- 
ally, by drawing triangles with shear sticks attached to 
the vertices, and flipping them). 



4. Parity-symmetric shear fields 

In this section we shall consider the consequences of 
parity-symmetry for the shear field, in particular with re- 
gards to the 3PCF. Consider first the example just men- 
tioned: take a pure E-mode field, e.g., coming from cos- 
mological ray-tracing simulation, and rotate all shears by 
45 degrees, to obtain a pure B-mode field. In the original 
E-mode field, peaks and valleys do not occur symmetri- 
cally: high-density peaks are present, but no deep valleys, 
since the dimensionless density contrast 6 is bounded from 
below by —1, but clusters of galaxies have a high density 
contrast (which, by the way, is the reason why clusters can 
be detected by shear measurements, but voids can not). 
Hence, the shear field will have tangentially-oriented pat- 
terns around mass concentrations, with no corresponding 
radial shear patterns present. The B-field, therefore, will 
have strong shear patterns with one circulation, but no 
corresponding ones with the opposite circulation. In other 
words, the B-mode field would carry circulation informa- 
tion, it would be possible to distinguish between a right- 
handed and a left-handed shear field, which obviously vi- 
olates parity-symmetry. Hence, the asymmetry between 
peaks and valleys in the convergence field k e from which 
the E-mode shear field is derived translates into an asym- 
metry between handedness of the B-mode field after the 
45 degree rotation. If the convergence field would be sym- 
metric, with as many peaks as valleys, or more precisely, 
if all odd moments of k would vanish (as is the case for a 
Gaussian random field), then this handedness problem in 
the B-mode field would not occur, and it would be parity- 
symmetric. What we shall show in this section is, that all 
correlation functions of the form 



( 7 ^(X 1 ) . . . 7 B (X m ) 7 f (X m+1 ) . . . 7r E (X„)) = , (12) 
if m is odd. 
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4.1. Parity transformation of E/B-mode shear fields Using these relations in (JTSJl, one obtains 



Consider a general shear field, which can be a combination 
of E- and B-modes. Assume that we flip the field along a 
line through the origin, which encloses an angle £ with 
the positive a^-axis. A point x is mapped through this 
flipping to the point x' = Ax, where the matrix A has 
the properties that det A = — 1, and AA = 1 is the unit 
matrix. This second property shows that 



x' = 4x; x = A x' . 
Explicitly, 

cos 2£ sin 2£ 



A = 



sin 2Q — cos 2Q 



(13) 



(14) 



The parity-reversed shear field 7'(x') = (P7)(x') can be 
derived as follows: Let the shear at position x be 7(x) = 
|7|e 2iv , where if is the angle the shear 'stick' encloses with 
the positive Xi-axis; the flipped shear is then oriented in 
the direction ip' which satisfies <p + p' = 2£, so that 



7 '(x') = | 7 |e 2i(2C - v) =e 4iC 7 *(x) 



(15) 



The same result can of course be directly obtained from 
JTJ, setting T = A. 

The shear field 7' will again be derivable from two 
fields ip E,Bl , for which the relation © is valid. Hence, 



2 \dxf dx'2 ) r v " ' dx' x dx' 2 



72(x') 



2 dx: 



The fields ip E - Bl will be related to the original ones by 

^ e '(x') = tt e ^ e (x) ; V B '(x')=7T B V B (x), (17) 

where the factors 7r E ' B = ±1 account for the possibil- 
ity that the ip are not scalar fields, but pseudo-scalars. 
Introducing ltT7|) into ((T5|) and using 

where summation over repeated indices is implied and use 
has been made of (| 1 31) . one finds 



7l (x') = — {A iX A jX - A l2 A j2 )^% - ir B A tl A^% , 
7 2 (x) = ~y (AnAji - A i2 Aj 2 ) ipjj + TT E AiiAj 2 'il> E ij (18) 

Using next the explicit representation l|14|l of the matrix 
A, one finds that 



{An A 



An A 



ip,u - ip, 



22 



cos4£ -I- ip 12 sin4£ ; 



AaA^.ij = - ^' 22 sin4C - V,i2 cos4C . (19) 



7 ((x') = 7T 




7T B ^ B 12 C0S4C 



sin 4C 



(20) 



= 71 (x) cos4( -(- 72 (x) sin4( , 

where in the last step we have used the transformation 
law (|15fl of the shear. Similarly, 

J f-J\ — /_B^ B l _ ^ B 2 _E„,,E 



7 2 (x') - tt* 



7r>* COS4C 



^ E ii-V^ E 



+ I 7T B ^ 2 +7T E 



22 



sin 4C 



(21) 



= 71 (x) sin 4( — 72 (x) cos 4£ . 

Comparing these last two equations with (Q, one finds 
that 7r E = +1, 7r B = — 1: hence, the deflection potential of 
the E-modc field is a scalar field, as expected, whereas the 
potential of the B-mode field is a pseudo-scalar, so that 

^ E '(x') = ^ E (x) ; V B '(x') = -^ B (x) , (22) 
or 

V/(x')=V*(x). (23) 

4.2. Consequences for parity-symmetric shear fields 

If a shear field is parity-symmetric, the two fields 7(x) 
and 7'(x) have the same statistical properties; in particu- 
lar, all their correlation functions are identical. If one has 
a pure B-mode shear field, this then implies that all odd 
correlation functions vanish: From 112211 one sees that the 
potential of the field 7' has the opposite sign of that of 7. 
Since the shear and the potential are linearly related, this 
implies that the signs of the two fields 7 and 7' are op- 
posite, rendering all odd correlation functions zero. More 
generally, this argument shows that if we decompose the 
shear field into an E- and B-field, as in @, all correlations 
of the form 

<7 B (X 1 )7 B (X 2 )7 B (X 3 )) = 0; 

<7 B (X 1 )7 E (X 2 ) 7 !(X 3 )) = . (24) 

Therefore, only the pure E-field 3PCF and the mixed cor- 
relator 

(7 B (X 1 )7 B (X 2 ) 7 !(X 3 )) 

are non-zero for parity-symmetric fields. A further gener- 
alization then yields the result l|12fl , namely that all corre- 
lators containing an odd number of B-mode shears must 
vanish. Furthermore, all odd correlation functions of k b 
and ijj B vanish, for the same reason. 

One example of this is well known from the 2PCF of 
the shear; the mixed correlator ( 7 t7x ) vanishes identically 
due to parity invariance; on the other hand, this corre- 
lation function is linear in ^k e k b ) (see Schneider et al. 
2002), and the result derived here applies. 
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4.3. Aperture measures 

Given that all third order statistics of the shear are lin- 
early related to the 3PCF, these results have further im- 
plications. Define the aperture measures for a point at the 
origin (Schneider 1996) as 

M ap (9) = /" d 2 tf Q(|tf|) 7 t(tf) ; 
Jo 

M y {0) = / d 2 tf Q(|tf|) 7x (tf) , (25) 
Jo 

where the tangential and cross components of the shear 
are taken with respect tot he direction of the center of the 
circle of radius 8, and Q(i9) is a weight function. These 
aperture measures have a number of very useful proper- 
ties. Employing the relation © between shear and the 
potentials in Fourier space, it is easy to see that 

M ap (0) = for pure B-field (26) 

and 

M x (0) = for pure E-field . (27) 

Therefore, these aperture measures cleanly separate E- 
and B-modes; they are therefore used to obtain this sep- 
aration in data fields, and essentially all detections of B- 
modes in second-order cosmic shear statistics have been 
made using (M 2 ) (van Waerbeke et al. 2001a, 2002; 
Hoekstra et al. 2002; Jarvis et al. 2003). Since Af ap is 
a scalar quantity, it can be conveniently used to define 
higher-order cosmic shear statistics; in Schneider et al. 
(1998), the third-order statistics was calculated us- 

ing quasi-linear perturbation theory of structure growth, 
whereas van Waerbeke et al. (2001b) calculated (M| p ) 
using the fitting formula of Scoccimarro & Couchman 
(2001) for the non-linear evolution of the bispectrum of 
the cosmic density fluctuations. Munshi & Coles (2003) 
derived this statistics using the hierarchical clustering 
ansatz, which is expected to give an accurate prescription 
of the cosmic density field on small scales. 

^Af^p) can be obtained from proper integration over 
the E-mode shear 3PCF similarly, (M^) can 

be obtained by integrating over (7^7,^7^ )■ The latter, 
however, vanishes identically for parity-symmetric shear 
fields, implying that (M^) = 0. For the same reason, 
(Af 2 p Af x ) = 0. For the aperture measures, the handed- 
ness argument used above becomes even more intuitive. 

Pen et al. (2002) employed third-order aperture statis- 
tics to the Virmos-Descart cosmic shear survey; they find 
a significant non-zero signal for all their four third-order 
moments. Their measured values of (M x ) and (Af 2 p Af x ) 
shows that the data is not statistically parity invariant. 
Hence, non-zero values of these quantitites cannot be ac- 
counted for by intrinsic galaxy alignments, as their correla- 
tion functions are expected to be parity-symmetric as well, 
nor to higher-order lensing effects - again, they produce a 
parity-symmetric shear field. The remaining explanations 
are that the cosmic variance is larger than estimated by 



Pen et al., or that there is a yet undetected systematics 
in the data. 

5. Discussion 

We have shown that for a shear field whose statistical 
properties are invariant under parity transformation, all 
odd correlation functions of the B-mode shear vanish iden- 
tically. This result has been derived by an explicit calcula- 
tion of the properties of the parity-reversed shear field, as 
well as by a simple consideration: Whereas an asymmetry 
between peaks and valleys in the convergence k e , which 
leads to non- vanishing odd correlation functions of the E- 
mode shear field, signifies that matter overdensities and 
underdensities behave differently, a similar asymmetry in 
the corresponding field k b yields an asymmetry between 
the two circularizations of the B-mode shear field, which 
for a parity-symmetric field is not permitted. In particu- 
lar, the 3PCF of the shear vanishes identically for a pure 
B-mode field. This result is at variance to some claims 
in the literature, and we have attempted to explain the 
reasons for this. 

One possibility to numerically generate a B-mode field 
which is parity-symmetric is the following: Take a pure 
E-mode field, as it is obtained from ray-tracing simula- 
tions, and rotate all shears by 45 degrees. Take a second 
realization of the same random field (i.e., a second ray- 
tracing simulation with identical parameters but different 
random number seed), rotate by 45 degrees and parity 
flip the shear, or, equivalently, rotate the E-mode field by 
—45 degrees. Then take the average of these two fields, 
which is then a pure B-mode field and parity symmetric 
(by construction). 

The result obtained here provides a very useful di- 
agnostics for potential systematics in cosmic shear data. 
Whereas real B-modes can be present in the data, as they 
can be generated by intrinsic galaxy alignments or higher- 
order lensing effects, e.g. coming from source clustering, 
all these effects should obey parity symmetry. A significant 
detection of a correlation function involving an odd num- 
ber of B-mode components, or a quantity derived from 
it (such as the aperture measures discussed above) is a 
sensitive probe of parity violation that probably can only 
be accounted for by systematics. It should be noted that 
a similar diagnostics has been used before for the 2PCF, 
namely checking that ( 7 t7x ) vanishes. 
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